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Abstract. We examine crossing probabilities and free energies for conformally invariant critical 
2-D systems in rectangular geometries, derived via conformal field theory and Stochastic Lowner 
Evolution methods. These quantities are shown to exhibit interesting modular behavior, although 
' the physical meaning of modular transformations in this context is not clear. We show that in 

many cases these functions are completely characterized by very simple transformation properties. 
In particular, Cardy's function for the percolation crossing probability (including the conformal 
' dimension 1/3), follows from a simple modular argument. A new type of "higher-order modular 

' form" arises and its properties are discussed briefly. 
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, 1. Introduction. There are extensive and well-known applications of modular invariance to 

On ' various quantities arising in conformal field theory (CFT), beginning with the work of Cardy 

[Cal] on the implications of modular invariance for the operator content of a given theory. How- 
I ever, in these cases the system is defined on a torus so that the modular invariance is implicit 

already from the definition of the problem. In this paper we examine the modular behavior of 
several quantities defined on compact two-dimensional figures with boundary, typically rectan- 
gles, where there is no obvious reason to expect modular properties. More explicitly, in the cases 
we consider, the modular operation S corresponding to r — > — 1/t is implied by a symmetry 
of the problem, but the operation T which takes r to r + 1 has no apparent physical inter- 
pretation; its origin is mysterious. Despite this, we are able display some new and interesting 
modular properties of known solutions to several problems from CFT and the recently developed 
^ I Stochastic Lowner Evolution (SLE) method [LSWl]. Examination of these properties for the 

' crossing probabilities (originally defined for percolation, but more recently extended to other 

^ processes via SLE) leads us to define a new type of modular quantity that we call an n*'* order 

modular form. Conversely, we show that by postulating a specific form for the crossing proba- 
bility (and similarly for the partition function), it is possible to reproduce the explicit forms for 
these quantities. 

In Sections 2 and 3 we briefly review percolation and the exact analytic forms of the crossing 
probabilities given by CFT [Ca2], [W] and transform these results into a form suitable for the 
present analysis. We also mention their generalization to other processes via SLE [LSWl]. 
Section 4 briefly recapitulates some of the basic properties of modular forms. In Section 5 we 
prove several theorems showing that the crossing probabilities are uniquely characterized by 
two very simple assumptions, a symmetry property and an assumption about the form of their 
g-expansion. In addition, we exhibit a possible generalization of one of the crossing probabilities 
to the SLE processes. In Section 6, we consider the (universal part of the) partition function for 
systems on a rectangle and show that the exact expression, already known by a CFT calculation 
[KIV], has a modular characterization of a similar kind. Section 7 considers the "higher order 
modular forms" that arise from the crossing probabilities. Roughly speaking, a first order 
modular form is an ordinary modular form and an nth order modular form is a function whose 
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deviation from modularity is an (n — l)st order modular form. Section 8 concludes the paper 
with a brief discussion. 

The subject we consider lies between mathematics and physics. Therefore some introductory 
matter from both fields has been included to make our treatment more comprehensible to those 
with background in one area but not the other. 

2. Crossing Probabilities. 

In this section, we first briefly review percolation (for a more complete treatment, see [Ke] or 
[StA]), and the definition of the crossing probabilities [LPPS]. Then we give the exact analytic 
forms for the latter quantities obtained via boundary CFT ([Ca2], [W]) and their generalization 
from SLE [LSWl]. The various crossing probabilities are displayed in three different ways: as 
an (ordinary or generalized) hyper geometric function of a certain cross-ratio A (to be defined); 
as the integral of an algebraic function of A; and directly in terms of the aspect ratio r of the 
rectangle (following [Zil] and [Zi2]). In the last case, the probabilities may also be written as 
series in rational powers of the parameter q = e^'^*'^, where r := ir. 

Percolation is perhaps the simplest non-trivial model in statistical mechanics. It is very 
easy to define, and exhibits a second-order phase transition between the percolating and non- 
percolating states. There are various types of percolation; for definiteness, we consider bond 
percolation on a lattice with square unit mesh. The system considered is a finite rectangular 
L X L' lattice where L,L' ^ oo with fixed aspect ratio r = width/height = L/L'. A bond is 
placed with (independent) probability p on each edge of the lattice. Consequently, there are 2^ 
possible bond configurations with < Nb < N, where Nb is the number of bonds in a given 
configuration and N is the total number of edges. The connected bonds in each configuration 
form clusters. Note that for any configuration, either there is a cluster spanning the lattice 
from left to right, in which case the dual lattice has no vertical spanning cluster, or there is no 
horizontal spanning cluster on the lattice and the dual lattice has a vertical spanning cluster. 
For p > Pc, as the lattice is taken to infinity, an infinite cluster appears. For p < Pc, there is no 
infinite cluster. In the case at hand, it is known that pc = 1/2. 

The quantities that we consider are universal, i.e., independent of the type of (isotropic) 
percolation and the lattice structure, as long as one remains at the percolation point Pc- In 
addition, they are believed to have conformal invariance properties. In particular, the crossing 
probabilities computed for two geometries which can be mapped onto one another by a conformal 
map should be the same. The universality and conformal invariance are not actually rigorously 
proven in the general case. However Smirnov [Sm] has recently proven that site percolation 
on the triangular lattice is conformally invariant in the scaling limit. This work also derives 
Cardy's formula and confirms the conjectured connection of SLE and percolation (see below). 
The universality (and conformal invariance for other percolation and lattice types) are supported 
by extensive numerical work and much other evidence. 

At Pc, the probability of a configuration including a cluster spanning the lattice from left 
to right is the horizontal crossing probability Uh. The probability of a configuration including 
a cluster connecting all four sides is the horizontal-vertical crossing probability Hhv These 
quantities depend only on the aspect ratio r because of conformal invariance. By the duality 
argument above, the horizontal probability, as a function of r, must satisfy 



n^(r) + n^(i/r) = i. 



(1) 



while by symmetry, the horizontal-vertical probability must satisfy 



Uh,{r) = n^l/r). 



(2) 



The conformal approach to the horizontal crossing proceeds by considering the Q-state Potts 
model, expressing the crossing as a difference of partition functions with certain non-uniform 
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boundary conditions, identifying the corresponding boundary CFT operators, and taking the 
hmit Q ^ I. Since the changes in boundary conditions are implemented by boundary operators, 
one finds that the crossing probabihty is given by a four-point boundary operator correlation 
function. The horizontal-vertical crossing is obtained by a similar procedure; the main difference 
is that the boundary operator is not the same as for the horizontal case. These derivations are 
fully described in the original calculations ( cf. [Ca2], [W] or [Kl]), and not particularly germane 
to our purpose here, so we omit further details. 

Instead of r, we can use an alternative parameter for the rectangles which is suggested by 
the conformal invariance property. If we choose a one-to-one conformal map from the rectangle 
onto the unit disk (or upper half-plane; it doesn't matter), then the only conformal invariant 
of the geometry is the cross-ratio A of the four points to which the corners of the rectangle are 
mapped. One then finds, by appropriate conformal manipulations, that the correlation function 
for horizontal crossing satisfies a Riemann differential equation with the two solutions F{X) = 1 
andF(A) = X^/^2Fi{l/3,2/5;4/3; X). One can pick the correct linear combination by imposing 
the physical constraints that F— >-OasA— >0(r— >oo) and F — > 1 as A — > 1 (r ^ 0). The 
result is Cardy 's formula [Ca2] 

^h{r) = A^/^ (1/3, 2/3; 4/3; A). (3) 

The cross-ratio A is given explicitly as a function of r in Eq. (12) below. The hypergeometric 
function appearing in Eq. (3) is very special, since its parameters a, b, c satisfy c — a = 1, so 
that it reduces simply to the integral of an algebraic function: 

This is a reflection of the fact that the hypergeometric differential equation satisfied by 11^ 
factors as 

>l-A))V3j.^ = 0. (5) 

For the horizontal- vertical probability lihvj the corresponding CFT analysis leads to the 
fifth-order differential equation 

for the function F defined by Ilh^(r) = F{X). The physical requirements of the problem are that 
^hv{r) be bounded, satisfy the symmetry condition of Eq. (2), which in terms of A translates 
into 

F{X) = F{1-X), 

and satisfy the asymptotic condition 

hm = , 

r^oo Uh{r) 

where 

Uh^{r) = Uhir)-Uh,{r) (7) 

is the probability of there being a horizontal but no vertical crossing. Applying these conditions 
to the differential equation (6) gives the explicit expression due to Watts [W] 

^hv{r) = ^ A 3i^2(l, 1, 4/3; 2, 5/3; A) , (8) 
^ 3 



where 3F2 is a generalized hypergeometric function. Again, because of the special form of the 
parameters, this has a simple expression as an integral: 

Un.ir) = f\u{l-u))-^'^dudt. (9) 

V OTT Jo Jo 

The differential equation (6) can also be written 

This form is of interest since 1, Hh, and Hhv span the solutions of the equation formed by letting 
the rightmost factor act on F alone, i.e. 

A(A(l-A))V3j.,A(l-A,)V3j.^ = 0. (11) 

Note that the full set of solutions of Eq.(6) (or (10)) is spanned by adding the functions log A 
and log(l — A) to the three just mentioned. Thus this problem is an example of logarithmic 
CFT ([Gur], [Fl], [Gab]). This behavior is to be expected on general grounds ([Ca3]), but has 
apparently not yet been explicitly exhibited for percolation crossing probabilities. However, 
[GL] calculates a closely related quantity. 

In what follows, it is convenient to consider the r-derivativcs Ii'f^{r) and 11^- (r) of Hh and 
Hhv- Note that these quantities are interpretable physically as probability densities: for instance, 
njj(r) dr is the probability that the rightmost point of any cluster attached to the left vertical 
side of an infinitely wide rectangle of unit height lies between r and r + dr [Zil]. Note that 
Eq. (11) reduces to second order when considered as a differential equation for the derivative. 

3. The Crossing Probabilities in Terms of the Aspect Ratio. 

In order to set the stage for an investigation of their modular properties, we next proceed to 
express 11'^ and 11'^- on the rectangle as functions of the aspect ratio r, using the classical result 
for the cross-ratio, namely A = A(ir) where A(r) is the classical modular function ("Haupt- 
modTiF') for the subgroup r(2) of PS'L(2,Z) . (All needed properties of modular functions and 
modular forms will be reviewed in Section 4.) This function can be given by many formulas, 

^ U ^ ..v(r/2)hj{2T)^ r?(r/2)i«r?(2r)« f Mr)^ 

where ??(r), t?2(''") and t?3(r) are the classical modular forms of weight 1/2 (Dedekind eta function 
and Jacobi theta functions) defined by 
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r,{r) = q'/^'llil-qn = ^(-l)"^(6"+i)V24 ^ ^1/24^1 _ ^ _ ^2 ^ ^5 ^ _ _ 
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^2(r) = ^g("+l)V2 = 2^§^ = 2q^/\l + q + q^ + ...), (13) 



nez 



Mr) = = , ,0 = l + 2q + 2q^ + 2q'' + ... 
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with q = e^'^^'^ and q := e"^'^ = ^/q. Note that A = 16q - 128^^ + 704^^ + . . . tends to like 
q and that A is a power series in g, not q. ( The appearance of q is typical for conformal field 
theory on the rectangle.) The derivative of A(r) is given by 

J-\'(r) = g ^(^/2)'N (2T)^^ 
2m ^ ' r]{T) 



_A'(r) = ^'JSL1^±J1^1>— = 8q - 128e + 1056?^ + . . . . (14) 



Now, since ^ = y^jyj-^^-, Eq. (11) can be rewritten 

d (A(l- A))V3 d (A(l-A))^/3 d 
Irr )J{7) 'dr A^(r) drr 

or, in view of Eq. (14), as the differential equation 

d r]{TY^ d 1 



F(A(r)) = 



dr 7?(r/2)8r?(2r)8 dr r]{TY 



fir) = 0, (15) 



for the function /(r) := ^F(A(r)). From this we can immediately write down two linearly 
independent solutions 

/l(r)=r/(r)^ 

(The factor — 27ri/3, of course, is just for convenience.) The function fi is a modular form of 
weight 2 (cf. Section 4). The function /2 can be decomposed as 

f2{T) = ^Mrf - Ww{r), (17) 

where i?2(t)'^, the fourth power of the theta function in (13), is an odd function of q and 

^ , , 1 16 2 364 3 13568 4 91614 . 

Mr) = -^,+ ^^<l' + ^^,' + ^^,' + j^^,' + - . (18) 

an even function of q. The function i?2 is again a modular form of weight 2, but fw is a new type 
of modular object whose transformation properties under the modular group will be discussed 

in Section 7. 

The r-derivatives of Ilh and Uhv can now be written in terms of /i(r) and /2(t) as 

97/3 2 

Kir) = -^1^;^^ hi^r) , UUr) = -8V3/2(zr) . (19) 
(See also [Kl].) The functions Hh and 11^^ themselves are then given by 

27/3^2 i-oo poo 

^hir) = ^^^^i I hiit) dt , Uh,{r) = SVs | /2(^^) dt . (20) 

Finally, we consider the recent generalization of Hh via SLE [LSWl] . This is a rigorous theory 
of stochastic conformal maps, driven by a Brownian process of speed k, B{Kt), which has been 
used to calculate the Brownian intersection exponents. For k > the hull of the process is 
generated by a path for alH > ([RS], [LSW2]). For < k < 4 the path is simple, while for 
K > 8 it is space filling. We will sec in Theorem 2 below that the limits k = 4 and k = 8 arise 
from modular considerations as well. The corresponding horizontal crossing probability is given 
by a generalization of Cardy's formula, 

FiX; .) = , ^^^:'/"^ A-V« (1 _ 1, 1; 2 - 1: A) . (21) 
^ ' ^ r(l - 4/Av)r(2 - 4/k) ^ kk k ' ^ ' 
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It is easy to show that F{X; k) satisfies the same duaUty condition (Eq. (1)) as Cardy's formula, 
and reduces to it when k = 6. ([Sm] proves that k. = 6 corresponds to percolation, as conjectured 
by [LSWl]). Further, the hyper geometric functions involved again satisfy c — a = 1, so that one 
has 

(This also makes clear where the normalizing constant in (21) comes from, since Eq. (22) and 
the beta integral give F(l; k) = 1.) 

There are apparently no SLE results for the horizontal-vertical crossing. Our theorem in 
Section 5 gives a candidate solution, at least up to one undetermined parameter, and includes 

the percolation case. 

By the same arguments as in the special case «; = 6, we can now write 

^FiXir),^) = ^ (r?(r/2)r?(2r))8-2^« ' ^'^^ 

where we have set a = 1 — 4/k for convenience. (The rhs of (23) has a ^-expansion beginning 
with a constant times q".) The modular properties of this function will be discussed in Section 5. 
Integrating (23), we find the formula 

n;,(r;a) = 24"7r^^ ^(^^)'°"''" (24) 

for the generalization Ilh{r;a) := F{X{ir),4/{l — a)) of Cardy's Ilh{r) = ^h{f", 1/3). 

4. Review of Modular Forms. 

Let Fi = (Z) be the group of 2 x 2 integral unimodular matrices, acting on the upper 
half-plane M = {r G C | ^(r) > 0} by r ^ 7(r) = ^ for 7 = ( ^ ^) G Fi. A modular form of 
weight A; G Z on Fi is a holomorphic function / : H — C which satisfies 

/(7(r)) = (cr + d)' fir) (r G H, 7 = C ^) € F^) (25) 

as well as a suitable growth condition at infinity (specifically, |/(r)| < C{y^ + y~^) for some 
constants C, ^ > 0, where y = 9(t)). A modular function on Fi is a meromorphic function 
satisfying (25) with A; = (i.e., simply invariant under the action of Fi on H) and a weaker growth 
condition at infinity, specified below. Every modular function can be written (in infinitely many 
ways) as a quotient of two modular forms of the same weight. One can also consider modular 
forms and functions on subgroups F C Fi of finite index, where (25) is required only for matrices 
7 G F. Other generalizations include allowing a character by including a factor v{'~f) on the rhs 
of Eq. (25), where |v(7)| = 1 for all 7 G F, or allowing A; to be a half-integer. In the case of 
the full modular group Fi , the collection of equations (25) can be replaced by the two equations 
/(r + 1) = /(r) and /(— 1/r) = r'^/(r), since Fi is generated by the two matrices T = (q J ) and 
S = "q"*^)- These matrices (or rather, the automorphisms of H which they represent) satisfy 
the relation 5^ = (^r)^ = 1. A second group which will play an important role for us is the 
theta group Tq consisting of matrices in Fi congruent to 1 or S modulo 2. It is generated by the 
two matrices S and T^. 

We can restate the modular invariance property (25) conveniently as f\kj = /, with the 
action "|fc7" of 7 = ( ^ ^) G Fi on functions on H defined by (/|fe7)(T) = (cr -|- d)~^/(7(T)). 
This shorthand notation will be convenient in what follows. 
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It follows from Eq. (25) that modular forms of a given weight are a vector space over C, 
denoted Mi^iTi) for forms defined over the full modular group Fi, or more generally Mk{T,v). 
The dimension of this space for simple groups F and small values of k is very small, a fact 
which leads to many non-trivial identities, a few of which we exploit below. Modular forms have 
power series expansions in non- negative powers of g' = e^'^*'^. Modular forms / over Fi whose 
g-expansion has no constant term, so that / — > as r — > oo, are called cusp forms. In general, 
cusp forms of weight k on subgroups F C Fi are defined by requiring that |/(r)| < CS(r)'^/^ 
for some C > and all r G H. For each weight k they form a subspace Sk{T) C Mfc(F). As 
examples, the spaces M2(Fti) and S'i2(Fi) are both one-dimensional, with generators i?3(r)'^ and 
?7(t)^'^, respectively, while A'(r) satisfies A'|2T^ = A', A'|2S' = —A' and hence is a cusp form of 
weight 2 on F^; with a non-trivial character. 

An important theorem for a modular form / G M/. (Fi , v) expresses a kind of "sum rule" on 
the total number of zeros. It may be obtained by integrating the logarithmic derivative of f 
around the boundary of a fundamental domain F = H/Fi (for more details see [Ko]). Let i^p{f ) 
denote the order of the zero of /(r) at the point P G H, while Vooif) denotes the exponent of 
the leading term in the g-expansion of /, i.e. Voo{f) = a if /(r) has a Fourier expansion of the 
form X^^o '^nff'^"^ with oq / 0. Then 

Note that i'i{f), i^pif), ^p{f) are non-negative integers, while fooif ) is integral if k is even and v 
is trivial but can be a rational or even real number in general. (If fooif) = then v{T) = e^'^*".) 
The points i and p = are fixed by the elements S and ST of order 2 and 3, respectively, 

and thus are differently weighted. A similar formula applies if / G Mk{T, v) for any subroup F of 
Fi, but the rhs is multiplied by the index of F in Fi and a different set of points and weightings 
appear on the Ihs. In particular, for / G Mk{Te,v), Eq. (26) becomes 

^oo(/) + ^i(/) + ^^.(/)+ E Mf) = 'l, (27) 

Peia/Te, PT^i 

where Vooif) is now defined by fooif) = a if / = Z^^o'^'^^"^" with ao / 0, since the local 
parameter at oo is g not g, and J^i(/) is similarly defined as the leading power of q in f\kT~^S. 
Again Vi{f) and vp{f) are (non-negative) integers, while i'i(/) and I'ooif) can be arbitrary 
(non-negative) real numbers. 

One can also look at meromorphic modular forms. Now vpif) both at "finite" points P G H 
and at cusps like P = oo or 1 are still required to be finite but may be negative. (For A; = this 
is the "weaker growth condition at infinity"' for modular functions mentioned at the beginning of 
the section.) The sum rules given above still hold in this more general context. The presence of 
poles is reflected in the growth of the coefficients of the g-expansion as follows: these coefficients 
have polynomial growth in the case of holomorphic modular forms, grow exponentially in the 
square root of the index if the function is holomorphic in H but has poles at the cusps (for 
instance, for the modular function A(t)), and have exponential growth if the function has poles 
at finite points. 

5. Modular Properties of the Crossing Probabilities. 

The purpose of this section is to show that the functions found by Cardy [Ca2] , Watts [W] , 
and Lawler et al [LSWl] as the solutions of crossing probability problems are characterized by 
certain very simple mathematical properties. This naturally raises the question (to which we do 
not know the answer) whether these properties can be seen by some a priori arguments for the 
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crossing probabilities, in which case our theorems would provide very simple new derivations of 
the results of these authors. Our results also include a possible generalization of Watts's formula 
for the "horizontal-vertical" crossing probability to the SLE processes. 

We will call a function 11 on the positive real axis a conformal block if it is expressible as a 
(real) power oi q = e~'^'^ times a power series in q. More precisely, if 



with a € R and ao 7^ 0, we call n(r) a conformal block of dimension a. This type of function 
is ubiquitous in CFT on a rectangle, appearing in partition functions and correlation functions 
as well. We will call n(r) an even conformal block if a„ = for n odd, so that n(r) equals 

times a power series in q = q^. Cardy's crossing probability lihir) satisfies this stronger 
condition. Notice that the convergence of the series in (28) for all r > implies that a„ = 0(c") 
for any c > 1 and hence the corresponding function in the upper half-plane (defined by the 
same series but with q replaced by e'^^'^) is holomorphic and is an eigenfunction of the operator 

: T I— > T -|- 2 or, in the case of an even conformal block, even of T : r 1— > r -|- 1. 

Note that the nomenclature "conformal" is only suggestive at this point, since this definition 
does not imply that n(r) is related to any CFT model. 

We are now in a position to state 

Theorem 1. Let n(r) be any function on the positive real axis such that 

(i) n(r) is an even conformal block with dimension a > 0; 

(ii) n(l/r) = 1 -n(r). 

Then a = 1/3 and n(r) is Cardy's function. 

Proof. Define {a„} by (28), and define -P(r) (for r G H) by the same expression as in (28) but 
with q interpreted as e"^*"^ rather than e~'^'^, so that P{ir) = n(r) for r > 0. It follows that 
P(r) is analytic in the whole upper half-plane. By property (ii) and analytic continuation, we 
have P(— 1/r) = 1 — P(r), while the fact that a„ = for n odd gives P(r -|- 1) = AP{t), where 
A = e'^*". Hence /(r) := P'{t) is holomorphic in H and satisfies f\2S = — / and f\2T = Af. 
Also, / is small at infinity because of the assumption a > 0. At this point there are two ways 
to complete the argument: 

A: From f\S = -/, f\T = Af we deduce f\{ST)^ = -A^f and hence, since (ST)^ = 1, that 
A^ = —1, i.e., a = m/3 where m is an odd integer. It follows that G 5'i2(ri). But 
this space has dimension 1 and is spanned by r/^^, as mentioned in Section A, so f = Crf 
for some C 7^ 0. The integration constant required follows from the condition (ii). 

B: From the modular properties mentioned, / G M2iT1.1V) for some character v. The rhs 
of Eq. (26) is therefore 1/6. Since each term on the Ihs is non-negative and fi, Up and 
Up are integers, the only possibility is Vi = Up = up = and 1^00 = a/2 = 1/6. It 
follows that A^ = 1 and that the quotient of by r/^"* is holomorphic, bounded, and 
Fi-invariant, and therefore constant. 

Consider the percolation crossing problem. It is interesting that if one assumes, following 
Theorem 1, that the horizontal crossing probability is given by a single conformal block, its 
evenness follows from the physics, since the boundary conditions used in the conformal analysis 
are the same on the two horizontal sides (see [Ca4] for further discussion of this point). 

Our next theorem generalizes Theorem 1 by dropping the assumptions of evenness and posi- 
tivity of a, but one has to add a growth condition on the coefficients of n(r) which was automatic 
in the even case. 



00 




(28) 
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Theorem 2. Let Hi (r) be any function on the positive real axis such that 

(i') ni(r) is a conformal block of dimension a G M with coefficients an of polynomial growth; 

(ii) ni(l/r) = l-ni(r). 

Then < a < 1/2 and Ili{r) = Ilh{r;a), the generalized Cardy's function of Eq. (24). 

Proof. We argue as in Method B above. First we define Pi{t) similarly to P{t). Note that 
it is a conformal block but not necessarily even. Its derivative /i(t) := Pi{t) is holomorphic 
and of polynomial growth because of the polynomial growth assumption on the a„, and satisfies 
/^|25 = -/i and /i|r2 = Afi, with A = e^^^^". Since S and generate Te, it follows that 
/i G M2{T0,v) for some v. We now apply Eq. (27). The rhs is 1/2. The assumption of 
polynomial growth implies that Vooif) and i'i{f) are non-negative real numbers, while fi(/) 
and vpif) arc non-negative integers. On the other hand, I'ooif) equals a if a 7^ and is > 1 if 
a = 0. Eq. (27) therefore implies that < a < ^, Voo{f) = ct, 1^1 (/) = ^ — a and all other vp 
vanish. Since the rhs of (23) has the same properties, these two functions must be proportional 
(their ratio is a function on a compact Riemann surface with no zeros or poles, hence constant); 
and then since ni(r) vanishes at infinity it must be proportional (and hence, by (ii), equal) to 
lift (r; a). 

We make a few comments about the interpretation of the lower and upper bounds a = and 
a = 1/2 in Theorem 2 from the modular and from the physics point of view. In the proof of the 
theorem, both bounds arose from the requirment that the modular form of weight 2 given by 
(23) be "holomorphic at the cusps," i.e., that the two numbers Vooif) and vi{f) in (27) should 
both be non-negative. If a < 0, then the leading power in the g-expansion of this function 
is negative and the integral in (24) diverges. If a > 1/2, then the integral converges and the 
function Xlh{r]a) still satisfies the functional equations in (i) and (ii), but it no longer satisfies 
the growth assumption: now vi{P[) < and the grow exponentially in ^/n as explained at the 
end of Section 4. (In Theorem 1 we did not have to explicitly make the assumption of polynomial 
growth because in this case there was only one cusp — compare equations (26) and (27) — so that 
the assumption a> already implied the holomorphy of the function /(r) = P'(r).) Prom the 
physics point of view the assumption a > is natural since the probability n(r) has to go to 
zero for large r, but the physical meaning of the polynomial growth condition is not obvious. 
However, both critical values have a physical meaning in terms of the SLE processes, as already 
mentioned in Section 3: the lower limit a = corresponds to the value k = 4 above which 
the hull of the process is no longer generated by a simple path, while the upper limit a = 1/2 
corresponds to the value k = 8 above which the path becomes space filling. At intermediate 
values, namely 4 < k < 8, the path is self-intersecting. 

Our final theorem in this section reproduces and generalizes Watts's formula for the percola- 
tion crossing probability 11^^. 

Theorem 3. Let a and Ili{r) be as in Theorem 2 and 112 (r) be a second function satisfying 

(iii) n2(r) = e~'^^^ J2^=o bne~^^^ for some /? G R, with of polynomial growth; 

(iv) n_(l/r) = n_(r), where n_ := Hi - 
Then 

(a) < /? < 1, /? 7^ a . 

(b) The function n_ (r) is given by the formula 

n_(r) = C{a,P) {n{it/2)n{2it))^-^^- J, (r/(m/2)r?(2m))8-24(/3-a) ' ^^^^ 



with 



(30) 



(c) // also (r) 0,'nd n_(r) are positive for all r > 0, then (3 > a and 



n2(r) = C(a,P) / , , \ / , , '} \ , dudt. (31) 

7r- (??(it/2)?7(2ii))8-24« 7^ (r/(m/2)?7(2m))8-24(/3-«) ^ ^ 

The functions Hi, 112 and n_ axe intended to be the generalizations of Hh, Ilhv and Hhv, 
respectively, and (31) agrees with the second formula in (20) in the case a = 1/3, P = 1. 

Proof. The argument again follows the proof of Theorem 1, Method B. First we define Pi(r) 
as in Theorem 2, and P2 and P_ analogously. Both these functions and their first derivatives 
/i, /2 and /_ are holomorphic and of polynomial growth in H. The modular transformation 
equations of /i were given in the proof of Theorem 2, while f-\2S = /_ and f2\T^ = -B/2, 
with B = e^'^*^. Thus the function v = /2//1 satisfies v\qS = 2 — v and v\T'^ = {B/A)v, so 
the function g := v' fi = /_(/{//i) — f'_ satisfies g\iS = g, and g\T'^ = Eg. But g is also 
holomorphic and of polynomial growth (because /i, the modular form given in eq. (23), has no 
zeros in H and only exponential growth at infinity), so ^ G M/^iVe^v) for some character v. We 
now apply Eq. (27) to g. The rhs is 1, so, since all terms on the Ihs are non- negative and all 
except = /3 and 1^1 are integral (the equation g\4^S = g implies that i'i{g) is even!), we must 
have 0</9<l, 1^1 = 1 — /3 and i'p{g) = for all P G H. The case (3 = can be excluded since 
then the expansion of the function P2{t) would begin with a constant term and its derivative 
/2 would have order (3+1/2 = 1/2 rather than (3 at infinity, and the orders of v and g at 
infinity would become 1/2 — a and 1/2 respectively, leading to a contradiction with (27). Hence 
< /3 < 1. The fact that g is modular of weight 4 on Vq and has the given orders of vanishing 
at all (finite and infinite) points now fixes it uniquely up to a constant: we must have 

9{t) = C, "^^^ 



(r?(r/2)r?(2r))i6-24/3 

for some non-zero complex number Ci (again, because the ratio of the functions on the left and 
the right is a meromorphic function with no zeros and poles and hence constant) and therefore 

= J= (r?(r/2)ry(2r))S-24(/5-a) (^2) 

for some non-zero complex number C2. To complete the argument, we must integrate Eq. (32), 
multiply the result by /i to obtain /2, integrate again, and then adjust the constants, if possible, 
so that all the conditions of the theorem are satisfied. 

We first note that the transformation equation v{t) + v{—1/t) = 2 implies that v{i) = 1, so 
that (32) integrates to 

Since also P'_{t) = /-(r) = /i(t) {1 — v{t)) and n_(r) is small at infinity (because both ni(r) 

and n2(r) arc), this implies formula (29) except for the determination of the constant C{a,(3). 
Next, we look at the behavior of the functions at infinity. The right-hand side of (32) has a 
g-expansion beginning C2 ^^""(1 (8 - 24/3 + 2Aa)q + ■■■). It follows immediately that (3, 
completing the proof of (a) of the theorem, because if a were equal to (3 then we would have 
v'{t) = C2 + 0{q) and hence v{t) = C2T + 0(1), contradicting the periodicity (T^-invariance) 
of V. If /9 > a, then v'{t) is exponentially small at infinity and we can integrate (32) to get 
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instead of (33). The constant here must be because the g-cxpansion of i; = /2//1 should 
begin with a multiple of g'^"", not . To get the value of the constant C2, we compare (34) 
(with C3 = 0) and (33), obtaining: 

' Ji (??(T72)r7(2T'))8-24(/3-«) 

Using Eqs. (12) and (14), we can rewrite this by setting u = A(r') as 

.1/2 



C2 / {u{l - u)f-''-^ du 
Jo 



-16"- 



But from the invariance of u{l — u) under u>-^ 1 — u and the standard beta integral we have 
This gives the formula 

= r(/? - a)2 

for 13 > a, and since the function defined by (33) must have a g-expansion of the form coq"~" + 
• g° + Ciq"~"~^^ + • • • for all values of /? (again, because of the requirement that v, and not 
merely v', be a single conformal block), it follows by analytic continuation that the same formula 
is true also for /? < a. We note, in any case, that under the assumption that n_(r) and 112 (r) 
are both non-negative for real r (which is certainly what we want in the physical situation, since 
these functions are meant to represent probabilities), we have < 112 (r) < and hence 

automatically P > a (let r ^ 00 !), as stated in (c). Finally, the value of the constant C{a,P) 
in (29) and (31) follows from the above formula for C2 together with the requirement that n_ 
and 112 should add up to Hi, whose normalizing constant was already given in Theorem 2, and 
analytic continuation to include the case P < a. This completes the proof of the theorem. 

The arguments in the last part of the proof say that the formulas (29)-(31) can be rewritten 
in terms of the variable A = A(r) as n_(r) = U^vi^', (X, P) and n2(r) = Ilhv{^;a, P), where 

n..(A;a,/3) = 2 ^^^^^^^ /'(t(l - t))""^ f'' {u{l - u)f—' du dt , 

U,,{X;a,P) = 2^^^^^^^ f\t(i-t)r-^ [\u{l - u)f-"-' dudt . 

These two functions add up to the function Hh{r;a) of (24) (or (22) with k, = 4/(1 — a)), and 
all three satisfy the differential equation 

J.(A(l-A))'+"-^A(,(i_,))i-.^f = 0, (35) 

generalizing Eq. (11). The space of solutions of this equation is spanned by 1, Tlh{r\a), and 
Iihv{r\a,P). 

Theorem 3 picks out for each a certain special functions n^^(r; a, /3) and Xlhv{T; a, /3), related 
to Tlh{r;a) in a nice way, by relatively simple axiomatic properties. The hope, of course, as 
already suggested by the notation, is that for a suitable value of /3 these functions really may 
give the correct horizontal-vertical crossing probability for the SLE process with k = 4/(1 — a), 
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although we have no real basis for this beyond its validity in the special case a = 1/3. (In 
particular, there seems to be no obvious reason from the physics why Yih^^j should be a single 
conformal block.) Even assuming that it is true, we are still left with the problem of determining 
or guessing what the correct value of /? is. For a = 1/3 we know that we must take /? = 1. 
If we make the simplest possible assumption that f3 is given by a linear function of a, and if 
there really is a /3 corresponding to a for every a in the interval (0, 1/2), then we must always 
have /3 = 1, because this is the only line through (1/3, 1) contained in the box (0, 1/2] x (0, 1] 
permitted by Theorem 3. Another special feature of the value /3 = 1 is that here, and for no 
other value of /? > 0, the differential equation (35) reduces to a purely hypergeometric one, so 
that we have the formula 

n^,(r; a, 1) = • A 3^2(1, 1, 2 - 2a; 2, 2 - a; A) , (36) 

TT 1 — a 

generalizing (8), with A = A(r). Just as in the discussion of the special case a = 1/3 in Section 2, 
the hypergeometric equation occurring here is highly degenerate and its three fundamental 
solutions have the simple form 

^ ^ {n-a)\ A"+" y> (n-2a)! A" 
' ^ n\ n + a ' (n — a)\ n 

n=l 



n=0 



(where x! := r(a; + 1)), with power series whose coefficients involve only two gamma- functions 
rather than six as would be the case for a generic 3F2 . 

Note that SLE corresponds to the Q-state Potts models via Q = 4 cos^ {Att / k) , for k > 4, 
as conjectured in [RS] (see also [BB] and [FW]). Thus quantities such as the horizontal- vertical 
crossing probability, specified in the Potts case, should also exist in SLE. 

6. Modular Properties of the Partition Function. 

In this section, which is directed more at readers conversant with the basics of conformal field 
theory, we briefly discuss a different situation in which modularity arguments can be used in 
statistical mechanics. The partition function of any conformally invariant system defined on an 
I X rectangle (with edges) is supposed to have a universal common factor Z{l,l'). Assuming 
the same conformally invariant boundary condition on each edge, it is possible to calculate this 
factor by use of CFT. Up to an undetermined multiplicative real constant C, the result [KIV] is 

Z{l,l') = Cl^/^r]{T)-^/\ (37) 

where r := /I) and c G M is the conformal central charge. The rhs of this is a modular form 
(of real weight, with character, and possibly with poles at infinity) on the full modular group 
Fi, so one can ask whether it is possible to reproduce Z from modular considerations. This we 
now proceed to do, using certain assumptions based on the physics of the problem, thus showing 
that modular properties of conformal quantities defined on rectangles appear in a wider range 
of problems than just crossing probabilities. The theorem below is a variant of an unpublished 
argument due to Cardy. 

First, it is obvious from the definition of Z that it satisfies (i) Z{l,l') = Z{l',l). Next, 
because of the "trace anomaly" effect due to the corners of the rectangle [CP], we know that Z 
is homogeneous of degree c/4, so we can write (ii) Z{l,l') = I'^/^dir) for some function (i(r). 
Finally, if we assume that only one conformal block contributes to Z, then we have (iii) d(r) = 
q-c/is QnQ^, for all I' /I > 0, where ao 7^ and q = e^'^*'^ as usual, because the leading 
behavior of Z as r ^ 00 follows simply from known results for the partition function on an 
infinite strip and because the fact that the boundary conditions are the same on all sides of the 
rectangle implies that any single conformal block must be even (by the arguments in [Ca4]). 
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Theorem 4. Any function Z{l,l') which satisfies the conditions (i), (ii) and (Hi) above must 
he given by eq. (37). 

Proof. The convergence of the series in (iii) for r G iM_|- imphes its convergence for aU r G H, so 
(i(r) is a holomorphic function in M. Substituting (i) into (ii) we find that the product /i(r) = 
r]{Ty/'^d{T) is invariant under r ^ — I/t (first for r/i real, and then by analytic continuation 
for all r G H), while (iii) implies that /t(r) is also invariant under r i— r + 1 and hence under 
the whole modular group Fi. Since assumption (iii) also implies that h is bounded at infinity, 
and since the group Fi has only one cusp, h is bounded on all of H/Fi and hence constant. 

Notice that if we had weakened assumption (iii) to just "(i(r) is a single even conformal block, 
of dimension a," then we could still deduce that a = — c/48 + 2n for some integer n > 0, and 
hence that Z is given by (37) if we assumed a < —c/48 + 2, while otherwise (i(r) would in 
general be the product of ?7(r)~'^/^ and a polynomial of degree < ra in the modular invariant 
function j{T). 

7. Eichler Integrals and Higher Order Modular Forms. 

In this section we describe in a little more detail the modular properties of the crossing 
probabilities studied in this paper. We will concentrate mostly on the functions Iih{r) and 

^hvii") of the original percolation problem as given in eq. (20), but analogous remarks would 
apply also to the functions occurring in Theorem 3 for values of a and /3 other than 1/3 and 1. 

In the theory of modular forms, the Eichler integral of a modular form /(r) of integral weight 
k > 1 is a {k — l)st primitive of /, i.e., a function /(r) in the upper half-plane whose (A; — l)st 
derivative is (a multiple of) /(r). If /(r) = Yl Q^~^" , then /(r) can be given explicitly by 
/(r) = [^{z — t)^~'^ f{z) dz or by /(r) = ^(n + aY~^an g""*"". These functions are no longer 
modular, but are "nearly modular" of weight 2 — k: if / is modular with respect to F, then 
(cr + d)'^~^ f {j{t)) is the sum of /(r) and a polynomial in r of degree A; — 2 for each matrix 
7=("^)Gr. In the special case k = 2, f{T) is (up to a constant) simply the integral of 

/ from r to oo and transforms via fi^ir)) = /(r) + C{^) for all 7 G F, where 7 1-^ (^(7) is 
a homomorphism from F to C. In certain cases (namely, when / is a Hccke eigenform with 
integral eigenvalues; we do not explain the details), the image of this homomorphism is a lattice 
A C C and the function / gives a map from the modular curve H/F to the elliptic curve C/A. 

In the case of the Cardy function, we see from equations (16) and (19) or (20) that n^(r) 
equals fi{ir), where /i(r) is the (suitably normalized) Eichler integral associated to the weight 2 
modular form /i(t) = ?7(t)^. The function /i gives a modular parametrization of the elliptic 
curve C/L, where L is the lattice spanned by 1 and ^ + ^VS- This elliptic curve has the 
Weierstrass equation Y'^ = + 1. The constant involving F(|)^ in front of the first integral in 
(20) is essentially the reciprocal of one of the "periods" associated to this elliptic curve. (Again, 
we omit details.) 

We now turn to the second function /2(t) in (16), which is a less familiar type of modular 
object: it is not modular, but its failure to be modular is given simply by multiples of the 
modular form /i . More precisely, we have 

/2(r + 2) = /2(r), r-V2(-l/T) = /2(r) - C/i(r), (38) 

where C = 2'^/^7r'^ /3T{l/3f . To see this, we observe that /2 is the product of /i with the Eichler 
integral /s of the modular form /3(r) = ?7(r/2)^?7(2r)^/r7(r)^^ of weight 2. The function f^ 
transforms under the generators of Fg by fs{T + 2) = oJ^f-iir) and r~^/3(— 1/r) = —fsir), where 
u; = e^'^*/^, so its integral /s transforms by /3(r + 2) = i^^/3(t) + ci and f^^—l/r) = — /3(t) + C2 
for some constants of integration ci and C2, the first of which is easily seen to be 0. Multiplying 
these equations by /i(r + 2) = w/i(t) and r~^/i(— l/r) = — /i(r), one finds equation (38). 
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Note that the transformation properties in (38) are the ones which were used and generaUzcd 
in Theorem 3 (in particular, the function occurring in (32) speciahzes to for a = 1/3, /3 = 1). 
Notice also that the second equation in (38), which can be written as f2\S = f2 — Cfi, can 
be combined with the transformation property = — /i to say that the linear combination 
/2 — |C/i is invariant under 5*. In terms of the original problem, the functions /i, /2 and 
/2 — are proportional to the derivatives of Ilh, U^^v and Hhv, respectively, and this last 

property is just a restatement of (the derivative of) equation (2). Finally, we can combine the 

two equations (38) by saying that the vector F = \ ^f] transforms under and S by 

^(r + 2) = (j;' j)F(r), F{-l/r) = (^Zc j)^^' (^9) 

i.e., it is a vector- valued modular form of weight 2 on the group Tq. 

Following this example, we define a second order modular form of weight k with respect to a 
subgroup r of SL2{1') to be a holomophic function /(r) which satisfies /|fc(7i — l)|fc(72 — 1) = 
for all 71,72 G F, or equivalently, if /|fe(7 — 1) is modular of weight k for all 7 G F (rather 
than being for all 7 as for an ordinary modular form). More generally, an n-th order modular 
form is a function satisfying /|fc(7i — 1) • • • |fc(7n — 1) = for all 71, . . . , 7„ € F or, in a fancier 
language, a function annhilatcd by the nth power of the augmentation ideal I = Kcr(Z[F] Z). 
We make a few general remarks about these higher order modular forms. (For further properties 
and examples we refer the reader to the recent preprint [CDO].) First, if we denote by M^"'\r) 
the vector space of nth order modular forms on F, then M^"^ (F) is always finite-dimensional 
and in fact of dimension at most (1 + r + • • • + r^~^)D, where D = dim Mfc(F) is the dimension 
of the space of ordinary modular forms of weight k on T and r is the cardinality of a set of 
generators of F. Indeed, if qi, . . . , G F are generators, then the map sending / to the r-tuple 
{f\k{ai - I),... ,f\k{ar - 1)) maps Mjf'\r) to (M^""^^(F))^ and has kernel Mfc(F), so the 
result follows by induction. This bound, however, is not sharp, since for instance for F = Fi 
(the full modular group, generated by the elements S and U = ST of order 2 and 3), the space 
M^"^(F) reduces to simply Mjt(F) for all n, the reason being that (1 — S*)" = 2"~^(1 — S) and 
(1 - C/)"(2 + [/)"-i = 3"-i(l - U) in the group ring Z[Fi], so that the equations /|(1 - 5)" = 
and /|(1 — U)"' = already imply f\{l — S) = and /|(1 — U) =0. The same argument 
applies to any group generated by elements of finite order. (In the situation studied here, the 
group involved is generated by elements of finite order, but the above argument can no longer 
be applied because of the presence of a character, i.e., because the diagonal terms in the two 
matrices in (39) are roots of unity but are not all equal to 1.) In general, however, the spaces 
Mj^\T) are larger than Mfc(F). In particular, one can construct non-trivial modular forms of 
(say) second order and weight 2 simply by multiplying a modular form of weight 2 by the Eichler 
integral of another modular form of weight 2, as was done for the function /2. Note, however, 
that /2 is an atypical second order modular form since for a general such function / we would 
only require that /|(1 — 7) is some modular form of weight 2 for each 7 G F, while in the case of 
/2 each of the functions /2|(7 — 1) is a multiple of the same modular form /i (i.e., the functions 
/i and /2 together are the components of a vector- valued but first order modular form F{t), as 
we saw above). 

In summary, the study of a problem coming from statistical mechanics has led to the consid- 
eration of a new and interesting type of modular object. 

8. Discussion. 

a. In this work, we examine the modular properties of crossing probabilities and their gen- 
eralizations, considered as functions of r = ir, with r complex. Consider, for definiteness, the 
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horizontal crossing probability n/i(r). The original problem is defined on a rectangle, which 
corresponds to a fundamental domain of the modular lattice Zr + Z. However, the function 
defined by analytic continuation of n^(r) is not the correct crossing probability for percolation 
on the parallelogram generated by 1 and t = ir when r is not real. This holds for all the crossing 
probabilities and their derivatives, and is a basic manifestation of the problem of understanding 
the connection between the conformal and modular properties at play here (unlike the situation 
for, e.g., the partition function defined on a torus, where it is quite natural to expect modular 
invariance). The physical quantity defined by Smirnov [Sm] is in fact given by (the real part of) 
^h{^) with r complex, but it also does not appear to have any natural modular properties. 

Another way to describe this issue is via the "mysterious" behavior of the crossing probabil- 
ities under the operation T, mentioned in the introduction. In our treatment this behavior is 
introduced by the assumption of one conformal block in Theorems 1-3. In CFT and SLE, on 
the other hand, it arises from the particular differential equations that the crossing probabilities 
satisfy. 

b. Note that 11^ for percolation (for which c = 0) satisfies the condition that the boundary 
conditions are the same on both horizontal sides (see [Ca2]) while Hhv does not, since the 
boundary conditions in that case are different on all four sides of the rectangle [W]. This explains 
why the former is an even conformal block (see [Ca4]), while the latter includes a function that 
is not. (In fact, for lih, two conformal blocks appear in the CFT calculation, but one of them 
is a constant, i.e. of dimension zero and with all coefficients vanishing except oq.) On the other 
hand, the block for 11^(1"', ct) for the SLE processes is not even except for a = 0, when it is a 
constant, or for a = 1/3, i.e. percolation, as a consequence of Theorem 1. It follows from the 
arguments in [Ca4] that this generalized crossing probability cannot be expressed as a difference 
of partition functions with boundary conditions the same on both horizontal sides unless the 
dimensions of the conformal blocks that appear differ by half-odd integers, since each partition 
function would necessarily be expressible as a sum of even conformal blocks. 

c. In some recent work, the plus-spin horizontal crossing function in the critical Ising model 
has been investigated numerically [LLS]. This quantity satisfies (1) and the simulations indicate 
that it is conformally invariant and vanishes as r ^ oo as e""'' with a ~ 1/6. Although its 
vanishing is consistent with the asymptotic behavior of n^(r; 1/6), further work [LSA] indicates 
that the agreement of 11^ (r; 1/6) with the numerical results is significantly worse than that of 
a particular Ising model CFT solution obtained by the authors. This solution is not given by a 
single conformal block, which is consistent with Theorem 2 (if polynomial growth is assumed). 
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